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ABSTRACT 

The collapse of baryons into extremely massive stars with masses > 10''^ Mq in a small fraction of 
protogalaxies at z > 10 is a promising candidate for the origin of supermassive black holes, some of 
which grow to a billion solar masses by z ~ 7. We determine the maximum masses such stars can 
attain by accreting primordial gas. We find that at relatively low accretion rates the strong ionizing 
radiation of these stars limits their masses to Af^ ~ 10'^ Mq (Macc/lO^"^ Mq yr~^)~, where Afacc is the 
rate at which the star gains mass. However, at the higher central infall rates usually found in numerical 
simulations of protogalactic collapse (> 0.1 Mq yr~^), the lifetime of the star instead limits its final mass 
to ~ 10^ Mq. Furthermore, for the spherical accretion rates at which the star can grow, its ionizing 
radiation is confined deep within the protogalaxy, so the evolution of the star is decoupled from that of 
its host galaxy. Lya emission from the surrounding H II region is trapped in these heavy accretion flows 
and likely reprocessed into strong Balmer series emission, which may be observable by the James Webb 
Space Telescope. This, along with strong He ii A 1640 and continuum emission, are likely to be the key 
observational signatures of the progenitors of supermassive black holes at high redshift. 

Subject headings: stars: formation - accretion - ISM: H II regions - cosmology: early universe - theory - 
galaxies: formation 



1. INTRODUCTION 

The existence of 10^ - 10^ Mq black holes (BH) in mas- 
sive g alaxies by 7, less than a billion years after the Big 
Bang ()Fan et al.ll2003l: IWillott et al.ll2003l : [Mortlock et all 
1201 lli . remains one of the great mysteries of cosmologi- 
cal structure formation. In the ACDM paradigm, early 
structure formation is hierarchical, with small dark matter 
halos at early epochs evolving into ever more massive ones 
by accretion and mergers through cosmic time. Hence, 
it is generally held that the supermassive black holes 
(SMBH) of the z - 7 Sloan Digital Sky Survey (SDSS) 
quasars grow from much smaller seeds at high redshifts. 
The origin of these seeds, and how they reach such large 
masses by such early times, remains to be understood. 
At least four main processes have been proposed for their 
formation (see IVolonteril[20T0t [Alexander fc Hickoxll2011t 
lKhlopov|[2QTol ): the collapse of Pop HI stars into 100- 300 
Mq BH at z ^ 20 (e.g . Madau & Rccs 20 0il:lAlvarez et all 
120091: lMilosavlievi6et al.ll2009t ITanaka fc HaimanI I2009D 
the direct collapse of extremely hydrogen molecule-poor 
primordial gas in ^ 10* Mq dark m atter halos int o 
1# - 10^ Mq BH at z ~ 10 (e.g. 'Bromm & Lqe 
20031 iKoushiappas etlll 120041: Begelman et al.. 200' 



shift and evolve into SMBH must account for how they 
become so large by z 7 and why their numbers at that 
redshift are so small, about 1 Gpc~'^. Pop HI seed BH 
are plentiful at z 20 - 30 but must grow at the Ed- 
dington limit without interruption to reach 10* - 10^ Mq 
by z 7. This is problematic because they and their 
progenitors either expel all the baryons from the shallow 
potential wells of the halos that crea t e them, so they are 
born starving" (IWhalen et al.l 120041: 1 Johnson fc Bromm 



2007HPelupessv et al.ll2007l: I Alvarez et al.ll2009l:lJeon et al 



Lodato fc Nataraian 
Omukai et al.l 120081: 



200a Spaans & Si M 12006 



20101 llnavoshi fc Omukail 12011 



Regan fc Haehnelt ,2009; ISethi et al.l 



al. 2003), 
into lO'' 



see also Colgate et 
the collapse of dense primeval star clusters 
10^ Mq BH (s ee e.g. iDiorgovski et al.l l2008t 



201 ih . or they eject themselves from their halos, and thus 
their fuel supply, at hundreds of km/s if they are born in 
core-collapse supernova explosions (Whalen & Fryer 2011 
in prep). Also, accretion onto Pop HI BH has been found 
to be ine fficient on small scales, typically at most 20% Ed - 
dington (iMilosavlievic et al.l l2009t iPark fc Ricottil |2011[ ). 
making the constant d uty cycle s required for sustained 
growth difficult (but seelLil l201lD . 

If halos can instead congregate into primitive galaxies 
of ~ 10* Mq at z ^ 10 - 15 that are devoid of the coolant 
molecular hydrogen (H2), they reach virial temperatures of 
~ 10'' K and begin t o atomically cool . Numerical simula- 
tions of this proces s (IWise et al.l 120081: iRegan fc Haehnell] 
120091: IShang et all I2OIOD find that in analogy to Pop 
HI sta r formation in much smaller halos at higher red- 
Umcmura 2001: Bromm et al. 2002 



iDevecchi fc Volonteril |2009() . and the collapse of primor- 
dial overdensi ties in the immediate aftermath of the Big 
Bang (see e.g. iMack et al.l [20071: iCarr et al.l[20Toh . 

The processes by which black holes form at high red- 



Abel et al. 


.2002: ,0'Shea fc Norman ,2007; .Yoshida et al. 


20081: Turk et al. 


120091: Stacv et al.ll2010l: Clark et al.ll2011: 


Greif et al 


2011 


), baryons rapidly pool at the center of 



the halo and form a hydrostatic object. But this ob- 
ject is thought to become far more massive than Pop 
HI stars modeled to date because atomic line cooling 
and the deeper potential well of the halo lead to much 
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higher infall rates at its center, 0.1 - 1 Mq yr~^ rather 
than the 10""* Mq yr~^ typical of primordial star-forming 
minihalos at z ~ 20. Such objects could collapse into 
blac k holes that are far mor e massive than Pop III BH 
(e.g. lShibata fc Shapiro|[2Q02l ). with Bondi-Hoyle accretion 
rates that allow them to grow into SMBH in less time (e.g. 
Wyithe & Loeb 2011; but see lDotan et~alll201lt) . 

One difficulty with this scenario is that it is not yet 
understood how primitive galaxies can form in the high 
Lyman-Werner UV backgrounds needed to fully sup- 
press H2 molecule formation and quench Pop III star 
formation in its constituent halos prior to assembly (1 
- 10 times tho se expected at the epoch of reioniza- 
tion at z ^ 6: iDiikstra et all 120081: IShang et al.] l2010t 
IWolcott-Green et al.ll20l'ir see also Agarwal et al. in prep) . 
Furthermore, it is not known if the object at the center of 
the protogalaxy even becomes a star, because at such high 
accretion rates it reaches very large masses on timescales 
that are short in comparison to Kelvin-Helmholtz times 
and the onset of nuclear burning (although Ohkubo et al. 
2009 have modeled the evolution of Pop III stars under 
accretion at much lower rates). It could be that the cen- 
tral object reaches such high entropies and densities that 
it is enveloped by an event horiz on before reaching the 
main sequence (|Frver et al.ll200ll ). If the central object 
becomes a star, what governs its rate of growth, its final 
mass, and thus the mass of the SMBH seed it becomes? 
Numerous authors have studied how stellar radiation regu- 
lates accretion onto Pop III protostars, but at much lower 
inflow rates th an those in the centers of collapsing pro- 



togalaxies (e.g. [S tahlcr et al. 1986; Omukai & Palla 2001; 



Qmukai fc Inutsuka.2002:.McKee fc Tan.2008; .Stacv et al., 



2011aHHosokawa et al.ll201ll) . Do the much higher accre- 
tion rates in primitive galaxies quench radiative feedback? 

We have performed extensive semianalytical calcula- 
tions of radiative feedback by Pop III supermassive stars 
in collapsing protogalaxies at z ~ 10. In § 2, we examine 
the primary forms of radiative feedback on accretion onto 
the star (deferring the processes that can be ignored to 
thorough examination in the Appendix). We also derive 
the maximum mass that the supermassive star, and hence 
the SMBH seed, can achieve as a function of accretion rate, 
accounting for its prodigious ionizing UV flux. In § 3, we 
estimate the maximum mass that the star can reach if it is 
limited by its finite lifetime rather than by radiative feed- 
back. We discuss the observational signatures of rapidly 
accreting supermassive primordial stars in § 4. Finally, in 
§ 5 we review the implications of final supermassive stellar 
mass for SMBH seed mass and the appearance of the first 
quasars in the Universe. 

2. RADIATIVE FEEDBACK-LIMITED ACCRETION 

We adopt an analytical approach to estimate the maxi- 
mum mass that an accreting primordial star can ultimately 
attain in a collapsing 10^ - 10* Mq protogalaxy. For sim- 
plicity, we assume that accretion onto the star is constant 
and spherically-symmetric. We first consider the ionizing 
UV radiation emitted by the star, which is the dominant 
form of radiative feedback on the infall and envelops the 
star with an H 11 region. In particular, we make the follow- 
ing assumptions on the nature of the radiative feedback on 
the accretion fiow (see Fig. 1): 



(1) Outside the H 11 region, gas is accelerated only by 
gravity and falls inwards unchecked by gas or radiation 
pressure. 

(2) Within the H 11 region, the gas is decelerated due to 
photoionization pressure only. This follows from the fact 
that the massive stars we consider here radiate at very 
nearly the Eddington limit. By definition, this implies 
that the force due to electron scattering in the H 11 region 
exactly cancels that due to gravity; therefore, for simplic- 
ity we can assume that the only net fo rce on the gas is due 
to ph otoionization pressure (see also lOmukai fc Inutsukal 
I2OO2D . 

In this Section we first determine how ionizing radiation 
from the supermassive star limits its growth via accretion. 
We then consider the effect of the Lya photons to which 
the majority of the energy in the ionizing radiation is con- 
verted. We examine and exclude other processes that can 
contribute to feedback on accretion in the Appendix. 

2.1. Suppression of Accretion by Photoionization 

As we are considering constant, spherically-symmetric 
accretion, the equation that links the number density n of 
hydrogen nuclei and the velocity v of the gas at a distance 
r from the star is: 

^^^•j _ Miicc 

47rr^u(r)/iTOH 

Here, we assume that the gas is primordial with an average 
atomic weight of fimn = 1.2mH, where toh is the mass of 
the hydrogen atom. 

Because the ionizing photons are trapped within the H 11 
region, we assume that the velocity of the gas at its edge, 
at a distance rnii, is just the free-fall velocity: 



v{rmi) = 



/2GMA 2 

V ^HII 



(2) 



We note that this velocity is much higher than the sound 
speed of the gas, at which infall is foun d to proceed from 
large radii in cosinological simulations ([Wise et al.l 120081 : 
IShang et al.]l2010l: iJohnson et al.ll201lD . so we neglect this 
small contribution to the gas velocity for simplicity. 

The radius of the H 11 region, rnii, is determined by the 
ionization front jump conditions, which require a balance 
between the rate Q at which ionizing photons are emit- 
ted from the star and the sum of the rate of hydrogen 
recombinations in the ionized volume (first term on the 
RHS in equation 3) and the rate at which neutral atoms 
enter the H 11 region through the accretion flow (second 
term on the RHS in equation 3). Following iRaiter et al.l 
(|2010D . we note that the rate at which ionizations occur 
is higher than Q by a factor P which accounts for cor- 
rections to the neutral hydrogen leve l populations at the 
high temperatures ( c^ 2-4 x lO'' K; e.g. lWhalen et al.ll2004l: 
I Alvarez et al.l 120061) expected for the H 11 regions of mas- 
sive Pop HI stars. To account for this rate, we take it that 
the total ionization rate in the H 11 region is Qcff = PQ, 
where P is the ratio of the average ionizing photon energy 
and the ionization potential for neutral hydrogen, 13.6 eV. 
Integrating over the density proflle of the gas in eq. 1 to 
find the radius at which the total number of recombina- 
tions in the enclosed volume balances the total number of 
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Fig. 1. — Schematic representation of accretion onto a supermassive Pop III star of mass Af» and radius r*. Gas accretes at a constant 
rate Mace, being acted upon only by gravity at r > thii, the radius of the H II region of the star. Inside rjjii the gas, whose infall velocity 
is v{r), is decelerated by momentum absorbed during ionizations, which is indicated by the shortening of the velocity vector arrows. While 
continuum photons with energies hu < 13.6 eV escape the H II region and do not couple strongly to the gas even outside thii, resonant line 
photons do couple strongly and so in principle can also impart momentum to the gas. In particular, a large amount of momentum is emitted 
by H I recombinations as Lyo; however, these photons couple so strongly to the gas that within r-Ly^ more than 50 percent of the emission is 
trapped in the accretion flow and cannot propagate outward and affect the dynamics of the flow. Within rjjii, Ly« photons are reprocessed 
into Balmer series photons which escape the H II region and may provide an observable signature of accreting supermassive stars (see §4). 



photoionizations, we arrive at the following expression for 
'"Hii: 

M 



ieS 



dr 



H 2 



47rr^w(r)^mH 



dr 



,(3) 



where ae 
coefficient 



1.4 X 

> 2 X 



10-" 
10'' 



cm 



^ is the recombination 
coemcient for > 2 x 10^ K photoionized primo rdial gas 
(jWhalen et al.' 2004': 'Osterbr ock &: Ferlandl l2006D . Follow- 
ing lOmukai fc Inutsuka ( 20021 ). for simplicity we neglect 
the second term in our calculations, as it is in general 
much smaller than Qcff- This is true in particular for the 
solutions that we find for the minimum accretion rate onto 
a star of a given mass; at much higher rates, this term may 
become large, likely limiting the radius of the H ii region. 

Finally, once the accreting gas crosses the stationary 
ionization front (as shown in Fig. 1), the equation that 
governs its deceleration due to momentum imparted by 
photoionizations is: 



dv 
dr 



a^hvn 



a^hv 



47rr^w(r)/xmH 



(4) 



where again we have assumed that the photoionization 
rate of hydrogen is balanced its recombination rate, aBn, 
and that the momentum transferred to the gas per pho- 
toionization is hiy/c. We note again that here we have not 



included the force due to gravity or that due to Thomson 
scattering, as these two forces cancel one another since the 
star is assumed to radiate at the Eddington rate. For the 
massive primordial stars that we consider here, we take 
the effective surface temperature of the star to always be 
~ 10^ K, which yields an average energy per ionizing pho- 
ton hiy ~ 29 eV. In turn, this implies P = hv/li.Q eV = 
2.1, which we assume throughout our study. In equation 
(4) we again choose case B, as we assume that recombina- 
tions to the ground level of hydrogen result in the emission 
of photons which do not deposit a net momentum to the 
gas (being emitted isotropically and contributing to the 
diffuse radiation field) Q 

To find the solutions for steady state accretion, we must 
solve the equation of motion equation 4. In particular, to 
find the minimum accretion rate for which an inflow so- 
lution exists, we search for solutions for which the infall 
velocity of the gas goes to zero at the stellar surface r*: 



w(n) = . 



(5) 



That is, we must simultaneously solve equations 3 and 4, 
under the constraints given by equations 2 and 5. These 
constraints yield the maximum mass that the star can 
achieve by steady accretion before its intense ionizing radi- 
ation halts infall at r > and prevents its further growth. 
On the other hand, they also yield the lower bound for the 
final mass of the central object, since for v{r^,) > radia- 
tive feedback fails to cut off accretion and the star will 
grow even faster than when gas merely comes to a halt on 



^ At the high densities in the H II regions we consider (Fig. 2) and with the enhanced population of neu tral hyd: 
temperatures in primordial H II regions l IRaiter et al.|[2010l) . a ~ 30 percent smaller 'extended case B' ||H 
coefficient may be more appropriate, but our final results are not strongly dependent on this choice. 



igen in the n = 2 state due high 
ummer &: StorevlllQSTI ) recombination 



4 



J. L. Johnson, D. J. Whalen, C. L. Fryer, and H. Li 




1 10 1 00 



Fig. 2. — The infall velocity (top panel) and density {bottom panel) of tlie gas accreting onto primordial stars as a function of distance r from 
the center of the star for three stellar masses: 10^ (dark blue), 10* (red), and 10'' Mq {light blue). Each star accretes at the minimum rate 
(as labeled) its strong ionizing flux permits, as given by equation 11. The dotted lines show the radius rjjn of the H II regions surrounding 
the stars, with their colors corresponding to those of the respective stellar masses. Outside the H II region, the gas is in free-fall, but after 
crossing into the H II region boundary it is rapidly photoionized and decelerated until it arrives at the stellar surface with velocity i;{r») = 0. 
The dashed lines show the trapping radius rLya for Lya photons from recombinations in the H II region (assuming a gas temperature of 4 X 
10"* K, see equation 22); partly because the vast majority of these photons are trapped within rLy^, we can neglect Ly« scattering feedback 
as discussed in § 2.2. Accretion at rates below those shown here is not possible, as in such cases run — > oo and photoionization pressure halts 
infall at all radii. 



its surface. 

The solution to equation 4 has the general form 



v{r) 



acc I — 1 

r — A 



(6) 



where is a constant whose value must satisfy the con- 
straint that ■(;(r*) = 0. This yields for K 



K 




(7) 



The equation for the infall velocity in the H II region then 
becomes 



v{r) 





(8) 

As shown in Figure 2, these solutions give a relatively con- 
stant velocity for the gas in the H II region until it is very 



close to the stellar surface. In turn, equation (1) implies 
that n{r) is essentially oc in the H ii region. 

Applying this v{r) (and hence n{r)) to the integral in 
equation 3 and evaluating it yields rjjn: 



3aB V 



(9) 



Finally, we must satisfy the constraint given by equa- 
tion 2 to find the solution for the minimum accretion rate. 
Substituting equations 8 and 9 into equation 2 and rear- 
ranging the terms gives the following quadratic equation 
in AP : 



3aB V ^ 



(10) 



whose solution is 



The Growth of Supermassive Stars and Black Holes 



5 



M2 

acc 



AG Ah. 



Noting that the second term in the square brackets is much 
larger than the first, we have 



Ma, 



3aB V 



(12) 



which exactly matches the accretion rate for which the H ii 
region breaks out to infinity (see equation 9) and therefore 
stops accretion entirely because photoionization pressure 
exerts a net outward force on gas at all radii Thus, we 
see that the H II region is confined and that accretion pro- 
ceeds at infall rates higher than those given by equation 
12 and that no infiow solutions exist for lower rates 

Now, by expressing Q ( where again Qeff = -PQ, with P 
= 2.1) and r* just in terms of the stellar mass M*, we can 
find the maximum stellar mass for which accretion onto 
the star is permitted (at Mace in equation 12). We can 
take it that, if the star i s thermally relax ecQ and radiating 
at the Eddington limit ()Begelmanll201Clt) . then the follow- 
ing two equations relate the mass M* of the star to its 
ionizing photon emission rate Q and radius : 



Q 1.6 X 10 



50 , 



100 M© 



(13) 



and 



= 3.7 R, 



■0 



M, 



He I photoionizations under the condition that within the 
H II region He i is also photoionized to He ii. For the 
relatively hard spectra of massive primordial stars this is 
a sound assumption because the number of He i-ionizing 
phot ons is coniparab le to the number of H i-ionizing pho- 
tons (|Schaereijl2002[ ). The ratio of the radiation pressures 
from helium and hydrogen photoionization is then 

PHcI "He ^'^Hel Q!b,HcI „n 

~ ~ 0.14 , (16) 

Pm Qb.hi 

where the first term is the ratio of helium and hydrogen 
number densities, which is ~ 0.1. The second term is the 
ratio of the average ionizing photon energies for He i and 
H I, which by integrating the stellar spectrum of a 10^ 
K primordial star we find to be 38 eV/29 eV ~ 1.3. The 
third term is the ratio of He i and H i case B recombinatio n 
coefficients, which is ~ 1.1 (jOsterbrock fc Ferlandll200"6t ). 
While this is a relatively modest increase in the radiation 
pressure, for completeness we include it as a multiplica- 
tive coefficient when solving equation 4. Other forms of 
radiation pressure are not so easily accommodated by our 
analytical approach but can be neglected without strongly 
affecting our conclusions, as we discuss in the next Section 
and in the Appendix. 

2.2. Trapping of Lya photons 

In the H II region of the accreting star, a large frac- 
tion of the ionizing radiation is reprocessed into Lyman 
a photons, which couple very strongly to primordial gas 
and could impart significant momentum to it, as they 
in princip le can scatter many time s across the H ii re- 
gion, (e g- 



Doroshkevichfc Kolesiiil] 119761: lOh fc HaimanI 



These are from Table 1 of Bromm et al.l ( 20011) and are 
in goo d agreement with ISchaereiT i 2002| ) and iBegelmanl 
(|20m ). With these expressions, we find that the solution 
to equation 12 is 



M* 



10^ Mp 



Ma, 







10- 



(14) l2002t lMcKe~ Tan 2001). However, while these authors 
showed that scattering Lya can strongly counteract the 
accretion fiow at relatively low accretion rates, at the ex- 
tremely high rates under which supermassive stars grow 
(e.g. Maee ^ 10~^ Mq yr~^) the impact on the accretion 
flow is reduced due both to the large momentum of the 
infalling gas and to the trapping of Lya photons within 
the H II region. 

To show this, we first compare the momentum in Lya 
photons to the momentum of the accretion flow. The for- 
mer IS given by dRaitcr et al. 2010[) 



^Lya 



PLya i^^A^ya 



lO^^gcms-2 



M* 



100 



where M^ is the maximum stellar mass attainable under 
accretion of gas at a rate Mace- 

We plot this maximum stellar mass in Figure 3. As we 
shall see, the feedback due to ionizing radiation emitted by 
the star can limit its mass at relatively low accretion rates, 
which we term 'feedback-limited accretion'. However, at 
higher accretion rates the flnite lifetime of the star governs 
its maximum mass; we term this 'time-limited accretion'. 
We discuss why this is so in § 3. 

Although we have focused on the pressure due to hy- 
drogen photoionizations, it is straightforward to include 

^ While in our solution the H II region formally extends to infinity in this case, for accretion to be terminated it need only extend to the Bondi 
radius, outside of which gas cannot be gravitationally captured by the star. 

^ Relatively small deviations from the solutions we have found for v{r) at rnii may result in a shock developing there. However, even for a 
strong shock for which there is a density jump of a factor of four, the minimum accretion rates in equation 11 change by < 20 percent. 

* At the minimum accretion rates we find for a given stellar mass (equation 15), the 'trapping' radius due to electron scattering (which is 
proportional to the accretion rate, e.g. IBegelmanl 119781 ) is alw ays smaller than the stellar radius given by equation 14. Therefore, the main 
sequence stars we consider here are thermally relaxed (see also lOhkubo et al.ll2009l '). At higher accretion rates, the radiation emitted by the 
star can become trapped due to electron scattering, resulting in an expansion of the stellar photosphere. In this case the star emits fewer 
ionizing photons than given by equation (13) and the radiative feedback is thus weakened, as we discuss briefly in the appendix. 

^ We note that in the high densities in the H II regions we consider (see Fig. 2), frequent collisions can prevent the radiative decay of hydrogen 
via two photon emission, which raises the Lya luminosity above that expected in the low density regime, where /3Lya — 0.68 fSpitzer,,197a) . 



where Llyq is the Lya luminosity due to recombinations 
in the H ii region, in which /3Lya (— 0-9) Lya photons, 
each with an energy hi^hya = 1-6 x 10~^^ erg, are emit- 
ted for every recombination (jOsterbrock fc Ferlandll20Q"6l : 
iRaiter et aLll2010[ )FI Here, we have also related the num- 
ber of ionizing photons emitted per second Q to the stellar 
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mass with equation 13. The momentum of the accre- 
tion flow is 

10 gems 



Ma 



6 X 



(l03 kms-i) 



M, 



where we have normahzed the infall velocity v to its typical 
value at the edge of the H ii region in the solutions shown 
in Fig. 2, since this is where Lyman a photons will be 
most strongly coupled to the gas (indeed, these photons 
are trapped within the accretion flow at thu, as shown 
below). Equating these two expressions and using the re- 
lation between Mace and M^^^i^^ in equation 15, we find 
the conditions under which Lyman a radiation pressure 
could halt accretion: 



< 600 km t 



Ma 



O.lMQyr- 



(19) 



Since the inflow velocity near the edge of the H ii region 
is well above this value at the typical minimum accretion 
rates we flnd (< 0.1 M© yr""'^; see Fig. 2), it is clear that 
Lyman a radiation pressure has only a small impact on 
accretion in comparison to photoionization pressure. In 
particular, over the range of stellar masses and critical ac- 
cretion rates that we consider (Figs. 2 and 3), we find 
that the infall velocity is ~ three times that in equation 
19. In other words, at the high accretion rates we find, the 
momentum of infiow is always roughly three times what 
could be countered by Lya scattering. 

It might be thought that if a Lyman a photon is 
scattered and then traverses the H ii region, it could 
subsequently be scattered many times, th ereby enhanc- 
ing the momentuETL it im parts the gas (see lAdamsl Il972t 
iMilosavlievic et al]|2009[ ). However, upon scattering from 
an atom in the rapid (~ 10'^ km s~^) accretion fiow, the 
photon would be strongly blue-shifted relative to the gas 
entering the opposite side of the H ii region; as a result, it 
would couple to gas much less strongly thereafter, greatly 
limiting the degree to which it could add more momentum 
to the gas (sec Dijkstra fc W yithc 201Q, on how galactic 
outflows enhance Lya escape by this process). 

Furthermore, another effect which dramatically lessens 
the impact of Lya photons on accretion is that their large 
resonant scattering cross section ensures that they are 
trapped in the H ii region and do not scatter o ut to larger 
radii. To see this, we follow the formula given bv lBegelmanI 
()1978[ ) for the trapping radius due to Thomson scattering, 
rtrap = MaccO'T/47rmHC. We define the trapping radius for 
Lya photons rj^ya by substituting the cross section for Ly- 
man a scattering CTLya for the Thomson cross section ctt 
for electron scattering, and by accounting for the fraction 
of neutral hydrogen atoms off of which Lya photons can 
scatter. This yields 



fhya 



5 X 



10^^ cm 



Ma„ 



/] 



HI 



104 K 



(20) 



X 10-14 cni2 (T/IO^K)-^ (e.g. IMilosavlievic et al.|[2009l ). 
The fraction of neutral hydrogen in the H ii region can be 
estimated by assuming photoionization equilibriurrlfl; we 
obtain 



/] 



HI 



QcffCTHI 



3 X 10-5 ( ^ ( ^ V 

VlOiOcm-sAlOiScm/ 



f Q. 



off 



V 1050 f 



(21) 



where cthi = 6 x 10^ cm is the photoionization cross 
section for H i (although it is likely slightly lower due 
to the relatively hard sp ectrum of the star; see e.g. 
iJohnson fc Khochfaij[20Tl[ ). Then with this /hi, Q from 
equation 13, n(r) from equation 1, and Mace in terms of 
M* from equation 15, we find: 



10 



12 . 



M, 



100 



-lO'^kms 



/ T 



VlO^K 



(22) 



Assuming a tempera ture of 4 x lO^ K for the photoion- 
ized primordial gas ( Whalen et al.l[2003 : iKitayama et al.l 
|2004|) , Lya trapping radii for three cases are shown in Fig. 
2. We see that rLya is ~ 2 - 6 r* and is smaller than rnii 
by roughly an order of magnitude. While this is a modest 
fraction of the volume of the H ii region, because of its 
steep density profile (n oc r"^) most of the Lya photons 
originate from this region. 

Integration of equation 3 shows that the fraction of re- 
combinations produced in the H ii region as a function of 



r is 



' HII 



jj). Consequently, the frac- 



,M0yr- 

where /hi is the neutral hydrogen fraction, T is the tem- 
perature of the ionized gas, and we have used (TLya = 5.9 

^ Comparing the rates of hydrogen photoionization and coUisional ionization using the rate given by ICenI II1992I) . we find that the former is 
much greater than the latter; thus, photoionization equilibrium is a vaUd assumption. 

^ The Boiidi rad ius is roughly rR„ ^ Hi ^ 10^'' (Af,/100 Mq) cm, assuming a sound speed in the protogalactic gas of 10 km s ^, following 



tion of Lya photons originating from within the trapping 
radius is (r^^ - r~^a)/{r-'^ - r^j^j), or > 0.5. Thus, most 
of them are trapped deep in the H ii region and cannot 
propagate outward to slow accretion at greater radii. Lya 
scattering thus removes only ^ 15% of the momentum of 
infall, not 30%, so we are justified in neglecting its impact 
on accretion flow. 

At this point, we can show that Lya photons created 
outside TLya but within rnii are also trapped in the accre- 
tion fiow at the boundary of the H II region due to the huge 
optical depth of the neutral gas to these photons. This is 
easily shown by equation 20, which for a largely neutral 
medium (i.e. /hi — 1) yields a trapping radius on the or- 
der of 10^1 cm for the lowest accretion rates given by our 
solution (equation 12); as this is much larger than both 
rnii and the Bondi radius, the maximum distance from 
which the protogalactic gas in can accrete onto the starQ 
Therefore, we can also safely conclude that Lya photons 
will not escape from the H ii region to affect the dynamics 
of the accretion fiow at larger radii. 

We note that although most Lya emission is trapped 
deep in the H II region, the thermal state of the gas may 
not be greatly affected because cooling vi a other atomic 
transitions in hydrogen still occurs (sec c.g'O mukai|[200ll: 
iRaiter et al. 2010; Schleicher et al. 2010). At the high den- 
sities within the H ii region atoms of both hydrogen and 



IWise et al.l pOOg) and lShang eTaLilmol '). 
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helium are readily excited to higher energy (n > 1) states 
by both collisions and absorption of photons; in turn, Ly- 
man series photons (principally Lya) are easily destroyed 
by ionizing these excited atoms (see e.g. Osterbrock & 
Ferland 2006). However, the radiative decay of excited hy- 
drogen atoms also results in significant emission of Balmer 
series photons as well as two-photon emission from the 
H(2s) state, to which the H ii is much less optically thick. 
The net result is that the energy in Lyman series photons 
is reprocessed largely into Balmer series and two-photon 
continuum emiss ion that escapes the H ii region and cools 
the gas (see e.g. ISchleicher et ani2010[ ). Thus, photoion- 
ization heating of the gas can still be balanced by efficient 
radiative cooling, and even though much of the stellar ra- 
diation is reprocessed in the H ii region it still eventually 
emerges and may be observable (see §4)0 It follows that 
the luminosity of the outgoing radiation is always the same 
as that of the star; since the star is taken to shine at the 
Eddington limit, this in turn implies that the force due to 
electron scattering indeed balances the force due to gravity 
everywhere in the highly ionized H ii region, as we assumed 
in our calculation. This also suggests that the tempera- 
ture of the gas is low enough to keep the flow supersonic 
so gas pressure cannot stop accretion, as we discuss in the 
Appendix. 

Overall, we conclude that it is the momentum imparted 
to the accreting gas by ionizing photons, despite their be- 
ing mostly converted to Lya photons that are in turn con- 
verted largely to Balmer series photons, that primarily reg- 
ulates the growth of supermassive primordial stars. 

3. TIME-LIMITED ACCRETION 

As noted by iBegelmanI ()2010t) . the lifetime of a very 
massive primordial star burning nuclear fuel at the Ed- 
dington rate ^Edd and growing by accretion at a constant 
rate is tyje ~ 4 Myr, twice the lifetime of non-accreting Pop 
HI stars of similar massif However, the stellar lifetime is 
in fact shorter than this for two reasons. First, nuclear 
fuel must be burned at a higher rate than the Eddington 
luminosity LEdd in order to support the star against its 
constantly increasing mass. Thus, the star consumes fuel 
at the rate given by the sum of the Eddington rate and 
the rate at which the binding energy of the star increases 
(Begelman 2010): 



in 



L 



GAUAL 



Edd - 



= L 



Edd 



1 



^trap \ 



(23) 



burning contributes to the gradual turnover of the stellar 
lifetime-limited maximum stellar mass curve shown in red 
in Figure 3. 

While equation (23) is valid when radiation can escape 
from the stellar surface at r* and the star is thus thermally 
relaxed, it ceases to be when radiation is trapped in the 
accretion flow within a radius rtrap above the star due to 
Thomson scattering. This occurs for accretion rates above 
Macctrap ^ 4.3 X IQ-^ Mg yr"! (M,/100 M0)3. As radi- 
ation is trapped in the accretion flow, the accreting mate- 
rial can not radiate its energy within rtiap. Upon reaching 
the stellar surface it therefore deposits energy at a rate 
G'M*Macc/''trap, instead of at the higher rate GAI^,A4acc/r* 
assumed in equation (23). Therefore, with this substitu- 
tion is equation (23), at accretion rates Mace > Mace, trap 
the star consumes fuel as twice the Eddington rate. This 
results in a decreased lifetime info, which in turn leads to 
the sharp turnover of the time-limited maximum stellar 
mass curve in Figure 3. 

The second reason that the stellar lifetime is shortened 
is also related to the trapping of radiation above the sur- 
face of the star, as whe n this occu rs the stellar envelope 
ceases to be convective ()Begelman| [2010')F"l This results 
in the envelope material never reaching the core and being 
available for fusion . The fuel supply of the star is thus lim- 
ited to the mass within its convecti ye core, whose fr action 
of the total mass is (equation 26 of lBegelmanll2010l ) 



Mc 



A'L 



0.54 



/ ''trap \ 



Ma 



f A'U 



Mq yr- 



V 106 Me 



(24) 

where we have assumed a central temperature of 10® K. 
This second effect comes into play on the right hand side 
of the dashed line separating the fully convective regime 
from the partially convective regime (also separating the 
regimes in which r^, > rtrap and r* < rtrap)- 

As shown in Figure 3, these two effects ultimately limit 
the final mass of the star to ~ 10^ at the highest 
accretion rates (~ 1 Mq yr~^) found in cosmological sim- 
ulations of the formation of supermassive stars via direct 
protogalactic collapse (e.g. Wise et al. 2008; Shang et 
al. 2010; Johnson et al. 2011). That said, we emphasize 
that our calculations have been done with significant sim- 
plifying assumptions, some of which we briclfy discuss in 
Section 5. In particular, we have adopted a simple model 
for the evolution of rapidly accreting supermassive stars; 
however, to be fully confident in our results will require 
dedicated stellar evolution calculations that account for 
the continued growth of the star over a large portion of its 
lifetime (Heger et al. in prep). 

Finally, we note that 'dark stars', which are pow- 
ered by dark matter annihilation, could grow to consid- 
erably higher masses than those powered by fusion (e.g. 
Freese et al.. 20 08). This is because the cooler surface tem- 



where r, is the radius of the star (equation 14), rtrap — 
Macca-T/47rmHC = 6 x 10^^ cm (Macc/M© yr~i) is the ra- 
dius at which radiation is trapped in the accretion flow by 
electron scattering, and Lsdd = 47rGM*mHc/crT = 1-2 x 
10"" erg s"i (M,/100 Mq) is the Eddington luminosity for 
a star of mass M* . This effect of super-Eddington nuclear 

* To be precise, for a star with an effective temperature of 10^ K, rougfily 60 percent of its luminosity is emitted in ionizing pliotons, which 
are later reprocessed into nebular continuum and recombination emission, such as Lya. 

^ This statement assumes that nuclear burning commences when the star has a mass much lower than its final mass. However, we note that 
for Pop III sta rs with final masses of ~ 100 Mq nuclear burning may not begin until the star acquires a substantial fraction of its final mass 
l|Omukai &: Pa lla 2003); if so, the lifetime of the star will be much closer to the ~ 2 Myr expected for > 100 Mq Pop III stars of constant mass. 

^'^ In this event, th e ste llar photosphere also expands to rtrap, as previous authors have found in the case of accreting primordial protostars 
l|Stahler et al.|[l986l : [Omukai & Palla 2003). Since the internal temperatures of primordial protostars are lower than on the main sequence, 
their opacities are somewhat higher than those due t o Tho mson scattering alone. Consequently, their radii are found to swell to more than 
rtrap and perhaps terminate accretion, as discussed in lOmukai fc Palla (.20031) . 
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Fig. 3. — The maximum mass to which a star can grow at a constant accretion rate Mace- Above the maximum mass Mt oc Mf/J given 
by equation 15 [blue line), radiative feedback shuts off accretion because the H II region of the star breaks out and prevents gas infall. For 
accretion rates > 10""*^ Mq yr"-*^, the maximum mass is set instead by the Hfetime of the star (red line). While at relatively low accretion 
rates the lifetime of an accreting massive primordial star is ^ 4 Myr, at high accretion rates the stellar lifetime can be shortened dramatically 
due to the super-Eddington luminosity at which nuclear fuel must be burned to support it during its growth. Also, at high accretion rates 
(right of the dashed curve), the star ceases to be fully convective and only the material in the convective core of the star is available for 
nuclear burning, further limiting the lifetime of the supermassive star. The dotted horizontal line marks the mass (~ 3 X 10^ M0) above 
which the final mass of a star is dictated by the limited time available for growth and below which it is governed by radiative feedback. The 
shaded area denotes the range of possible stellar masses for accretion at a constant rate. 



peratures of such stars exert much less radiative feedback 
on accretion and dark matter fuel may last for a much 
longer time than nuclear fuel. However, recent high reso- 
lution cosmological simulations show that primordial stars 
forming at the centers of dark matter halos may not remain 
there because of dynamical interactions, as required for the 
continual capture and annihilation of dark ma tter in their 
interiors (jStacv et al.ll2011bHGreif et al.ll201l[ ). Therefore, 
it may be that the final masses of dark stars are not so dif- 
ferent fr om those expected for standard primordial stars 
(see also iRipamonti et al.ll20l"ot) . If, however, dark stars 
do grow to be more massive they would exhibit spectral 
properties distinct from t hose of solely fusion-powered su- 
perm assive Pop III stars (IZackrisson et al.ll2010t lllie et all 
l201lh . as we discuss next. 

4. OBSERVATIONAL SIGNATURES 

As supermassive stars are intense sources of radiation 
that could be detected by current and future surveys, we 
now examine their observable signatures. The fact that ra- 
diative feedback in most cases cannot terminate accretion 
onto supermassive Pop III stars in collapsing protogalax- 
ies implies that their H II regions will be confined deep in 
their host halos for most of their lives. It follows that their 
ionizing radiation is largely reprocessed into nebular emis- 
sion instead of escaping the halo and rcionizing the IGM. 
As discussed in § 2.2, because the accretion flow is opti- 
cally thick to the Lya photons, they cannot directly exit 
the halo and be observed. These photons are instead fur- 
ther reprocessed largely into Balmer series photons which 
do escape the H II region and halo and propagate into the 



IGM. Therefore, one likely signature of rapidly accreting, 
isolated supermassive Pop III stars in high-redshift pro- 
togalaxies is a strong Balmer line flux accompanied by a 
conspicuous lack of Lya emission. A detailed radiative 
transfer calculation is necessary to quantitatively predict 
luminosities for all the Balmer lines in hydrogen, but we 
can place lower limits on the Ha flux, which is expected 

to be the do minant li ne. 

Following [Schaererl (|2002l) and iRaiter et al.l (|2010[ ). we 
compute the luminosity Lho in Ha as a function of QeS as 
discussed in §2.1. Relating this ionizing photon emission 
rate to stellar mass with equation 13 yields Lho — 4 x 
10'^* (Af*/100 Mq) erg s~^; we again emphasize that this 
is a lower limit because we exclude reprocessing of Lya 
into Ha, which may dramatically boost Ha luminosities 
above those estimated here. The Ha flux at redshift z is 

jHa — 



10" 



erg s ^ cm ^ 



1040 erg s-i 



1 + z 
10 



(25) 



where Di,{z) is the luminosity distance to redshift z. We 
plot this flux as a function of stellar mass and redshift 
in Figure 4, with the arrows on the Ha curves signifying 
lower limits. 

The hard spectrum of ionizing radiation from hot pri- 
mordial stars also creates an He ill region from which 
a large luminosi t y in He 11 A 1640 emiss io n is expected 



(iQh et all [2001 



_ _ . iTumlinson et al. 

I2OOII: ISchaereij|2002t iJohnson et al.ll2009[ ). Since the opti- 
cal depth to this line is low, we also estimate how much of 



2OOII: iBromm et all 
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Fig. 4. — The flux in Ha (top panel) and He II A1640 (bottom panel) from accreting main sequence primordial stars, as a function of 
redshift z, for stellar masses of 10^ (dark blue), 10^ (fed), and lO'' Mq (light blue). The dotted black lines show flux limits for Scr detection 
of these lines in 10* second exposures with the NIRSpec and MIRI instruments on the JWST, operating at resolutions R = 1000 and R = 
1200-2400, respectively. NIRSpec operates at wavelengths of ~ 1 to 5 fim, while MIRI covers the range from ~ 5 to 28 fim. As discussed in 
§4, as we have not accounted for the reprocessing of trapped Lya photons into Balmer series photons (including Ha), the Ha fluxes shown 
here are lower limits, as indicated by the arrows. While detection of Ha from stars with masses of > 10* Mq may be possible out to very 
high redshift (e.g. z > 10), only accreting primordial stars with masses of at least 10^ Mq are likely to be detectable in He II A1640. 



its flux exits the halo. For the large ionization rates and 
densities in our stu dy, we can apply the standard model of 
iRaiter et al.l (|2010l) to compute the luminosity of this line, 
which we find to be Li64o ^ 2 x lO^^ (M^ /100 Mc^) erg 
s~^, w here we again have used Table 1 of iB romm et al.l 
(|2001h for the number of He ii-ionizing photons produced 
as a function of stellar mass. Then, replacing Lho with 
Li64o in equation 34, we derive the flux in the He ii A 1640 
line, which we plot with fna in Figure 4 as a function of 
stellar mass and redshift. 

Also shown in Figure 4 are detection limits for two in- 
struments that will be on board the James Webb Space 
Telescope {JWST): the Near-Infrared Spectrograph (NIR- 
Spec) and the Mid Infrared Instrument (MIRI) . The black 
dotted curves in Figure 4 denote the flux limits for 3<t de- 
tection of the Ha and He ii A 1640 lines in an exposure 
time of 10^ s at resolutions R — 1000 for NIRSpec and R 
= 1200-2400 for MIRI (Gardner et al. 2006) E3 Because 
ionizing photon rates vary linearly with stellar mass, so 



do recombination line fluxes. Consequently, only the more 
massive accreting primordial stars will be detectable by 
the JWST. Noting again that our Ha line fluxes are lower 
limits, stars with masses > 10^ Mq may be detectable in 
Ha out to z > 10, while accreting stars with masses on the 
order of 10^ M0 could perhaps be detected out to some- 
what lower redshift. Detection of the He II A1640 line will 
only be possible from stars with masses of at least 10^ Mq 
for sufficiently long exposure times. 

We note that the enhanced Balmcr line emission from 
these objects may pose some difficulty to their identifica- 
tion as supermassive Pop III stars. As mentioned earlier, 
the large ratio of He 11 A1640 fiux to hydrogen recombina- 
tion line fiux in Lya or Ha is unique to Pop HI stars. If, 
however, much of the Lya emission is converted to Ha, the 
drop in this ratio could mask the primordial nature of the 
star. That said, the absence of any Lya emission would 
still likely reveal the central object to be an accreting su- 
permassive star. Another attribute of rapidly accreting su- 



JWST bandpass data can be found at www.stsci.edu/jwst/science/sensitivity. 
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permassive stars is strong continuum emission below the 
Lyman hmit, which is the sum of the stellar conti nuum 
and the nebular continuum (e.g. iRaiter et al.l l20ldl ): the 
latter would be substantial due to the exceptionally high 
densities expected in the H ii regions of these stars. In- 
deed, it is possible that the sum of the stellar and nebular 
emission could be detected, for example, in the Deep- Wide 
Survey to be carried out with the Near-Infrared Camera 
(NIRCam) on the JWST. 

An important obstacle to finding these objects is that 
their numbers could be small. As discussed by previ- 
ous authors (|Bromm fc Loebl l2003t iDiikstra et all |2004 
120081 see also Agarwal et al. in prep), supermassive 
stars cannot be so abundant that the black holes they 
create exceed observed limits on the black hole mass 
density and x-ray background. This suggests that find- 
ing the black holes may be easier than detecting their 
progenitors. First, they could accrete material for at 
least 10^ yr and be luminous for far longer than the 
stars that created them. Second, future x-ray missions 
such as the J oint A strophysics Nascent Universe Satel- 
lite (J ANUS) (iRoming 2008: Burrows et al.. 20101 L6>^- 
STER (|Gorensteinll201lh . SVOM (jCotz et al.Tl200l . and 
the Energetic X -ray Imaging Survey Telescope (EXIST) 
(jGrindlavl |2005[ ) will perform all sky surveys with far 
greater coverage than the JWST. However, if the nascent 
black hole does not have an accretion disk i t may only emit 
weak l y in x-rays at birth ( Frver et 311120011 : iFrver &: Hegeij 
[ml [s uwa fc lokal 1201 ll but see also Komissarov fc 
Barkov 2010). If so, the SMBH seed does not become visi- 
ble until it begins to acc r ete surrounding protoga l actic gas 
(iHaehnelt fc Reej|l993t iKuhl en fc Madaul l2005t [Li et all 
I2007t IVolonteri fc Begelmanl l2010: Jo hnson et atlboilD . 

Supermassive stars could also be detected if they 
explode as luminous supernovae. However, previ- 
ous studies have concluded that > 10^ Mq Pop HI 
stars collapse to black holes without an explosion 
(iFuller et al.l 119861: iFrver fc Hegeil 120111; iMontero et all 
I2011h (but see lOhkubo et al.ll2006D . At 14 - 260 Up, 
howe ver, pair-instability supernovae occur (jHeger et al.l 
I200a) and may be observable by fu t ure missions such 
as the JWST (IWise fc Abe l' '200 ^ IScannapieco et all 
20051: IWeinmann fc Lillvl 12005: Jogge rst fc WhalenI 120111 



Kasen et al. ||201lD. Finally, we note that much of the con 



tinuum and line emission from rapidly accreting super- 
massive stars will appear in the near infrared background 
(NIRB) today. Although their contribution to the NIRB 
may be small if they are rare, it could be detected by mis- 
sions such as the Cosmic Infrar ed Background Experiment 
(GIBER; e.g. iBock et al.ll2006[ ). which is designed to find 
signatures of primordial galaxy formation at z > 10. 

5. DISCUSSION AND CONCLUSIONS 

We find that the masses of the stellar seeds of SMBH 
forming from baryon collapse in early protogalaxies are 
primarily governed by their intense ionizing UV flux and 
their finite lifetimes. For spherically-symmetric accretion 
at constant rates Mace ^ 0.1 yr"-'^, the maximum mass 
the star can reach is governed by radiative feedback and is 
~ 10^ (A/acc/10"^ Mq yr-i)? Mq. At higher masses, 
the H II region breaks out to large radii and terminates 
accretion. We have verified that other forms of feedback. 



such as gas pressure, radiation pressure from trapped line 
emission, photodissociation of H~ , and scattering of Lya 
photons are much less effective at slowing accretion (see 
§ 2.2 and the Appendix). 

At accretion rates above > 0.1 Mq yr~^ the lifetime of 
the star, not radiative feedback, determines its final mass 
by limiting the time for which gas can accumulate on the 
star. Radiative feedback limits supermassive Pop HI stars 
to final masses of ~ 3 x 10^ Mq and time-limited accre- 
tion limits them to ^ 10^ Mq at the highest accretion 
rates (~ 1 Mq yr~^) found in numerical simulation s of 
protogalactic collaps e (jWise et a l. 2008: Shang ct al. 20101: 
I Johnson et aIll201lD . 

We caution that our analytical calculations do not 
account for all conceivable effects that c ould stem the 
growth of supermassive stars. For instance. iMcKee fc TanI 
(2008) note that rotation of infalling gas leads to lower 
circumstellar densities and large r H ii regions with 
greater radiative fe edback (see also iHosokawa et 311120111 : 
IStacv et"ari l2011al) . implying lower final stellar masses. 
Our spherically-symmetric calculation excludes rotation, 
so the mass limits we find for a given accretion rate 
(Fig. 3) ar e upper limits. We do, h owever, find agree- 
ment with (jOmukai fc Inutsukal |2003) . who performed a 
similar calculation, although they only considered the 
growth of stars up to ^ 10^ Mq. Furthermore, we also 
note that the spherical symmetry and constant accretion 
in our models do not address accretion that is episodic 
or clum py and self-shielded from ionizing radiation from 
the star (IWhalen et al.ll2008l [20T0I : iKrumholz et al.ll2009D . 
These processes could cause accretion to proceed at lower 
time-averaged rate s than those predicted here (but see 
iKuiper et ani20TTI) . 

Another process that could truncate the growth of su- 
permassive stars well before feedback or stellar lifetimes 
is the onset of a general relativistic instability in the 
core of the star that caus es it to collapse when it be- 
comes sufficiently massive (lGhandrasekhad[T96l . Such 
instabilities are predict ed to set in once the sta r has 
grown to ~ 10^ Mq (jlbenl 119631: iFowleij I1964D . but 
stellar rotation coul d stabilize it against collapse up to 
much larger masses (lFowleij[l966l : [Bisnovatvi-Koganlll998l : 
iBaumgarte fc Shapiro! 119991 )^ While our models ignore 
rota tion, the accreting gas is likely to have some (e.g. 
iGolgate et al.ll2003i ): if so, the supermassive star will in- 
herit the angular momentum of the gas from which it 
formed and perhaps bypass the general relativistic insta- 
bility. Additional insight into the processes that limit the 
growth of supermassive stars will be gleaned from stellar 
evolution calculations accounting for the continual accre- 
tion of mass at high rates (Heger et al. in prep). 

If supermassive stars formed and grew to masses of > 
10^ Mq in the early universe, recombination emission from 
their H 11 regions may be bright enough to be detected 
by future missions such as the JWST. In particular, Lya 
photons trapped in the accretion flow are reprocessed into 
B aimer series photons that could escape into the IGM. 
Gonsequently, the formation of these objects is accom- 
panied by distinctive strong Ha emission together with 
strong continuum and He 11 A1640 emission, the latter aris- 
ing from the hard spectrum of hot Pop HI stars. However, 
their black holes may be easier to discover in observational 



The Growth of Supermassive Stars and Black Holes 



11 



surveys, given the small numbers and brief lifetimes of 
their progenitors. Indeed, these black holes may be the 
very ones that have already been found at the centers of 
massive galaxies and quasars at high redshifts. 

Current numerical simulations of SMBH seed formation 
in early protogalaxies are now at an impasse because the 
formation of the hydrostatic supermassive protostar re- 
stricts Courant times to values that are too short to evolve 
central flows for even one dynamical time (but see Johnson 
ct al. 2011 for an alternative approach using the sink par- 
ticle technique). Our results suggest that in the early and 
intermediate stages of the growth of the star its evolution 
is essentially decoupled from flows on even slightly larger 
spatial scales. Consequently, it should now be possible to 
retreat from the extreme spatial resolutions previously ap- 
plied to the protostar and evolve flows at the center of the 
galaxy over enough dynamical times to capture its struc- 
ture at the time of the death of the star and breakout of 
x-rays from the BH seed into the IGM. Thus, it will soon 



be possible to witness the births of the first quasars in the 
universe with supercomputers. 
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APPENDIX 
NEGLECTED FEEDBACK PROCESSES 

Here, we discuss several processes which have a neghgible impact on our estimate of the the final maximum stellar mass 
given by equation 15. 

He II Photoionization 

To assess the relative importance of He ii photoionization in slowing the accretion flow in the He ill region, we follow 
the argument in Section 2.1 for He i photoionization. For the ratio of He ii photoionization pressure to H i photoionization 
pressure, we have 

PHcII niic hvHcII "B.HcII -. , ..,x 

— r — ^--"^ I [^^) 

Pm "•H hmi ttB.Hi 

wher e we have taken it that hvncw / hv-ni — 54.4 eV / 29 eV ~ 1.8, and that aB.Hcii/aB.Hi — 6.4 (jOsterbrock fc Ferlandl 
120061 ). Therefore, within the He III region, the radiation pressures due to He ii and H i photoionizations are comparable. 
This is convenient because in t he He III region it is recomb ination emission from He ii, not stellar photons, which largely 
keeps hydrogen photoionized (jOsterbrock fc Ferlandl [20061) . Therefore, while this diffuse recombination emission does 
not transfer appreciable outward momentum to the gas, the stellar photons that ionize He ii transfer roughly the same 
momentum to the gas that they would in ionizing Hi. As a consequence, solving in detail for the pressure due to He il 
photoionization would result in an almost identical solution to the one we have found by just treating the H ii region. 

Photodetachment 

In order for very high accretion rates onto a super massive star in a primordial protogalaxy to be realized, the H2 fraction 
in the accreting protogalactic gas mu st be very low (jBromm fc Loebll2003t ITodato fc Nataraianll2006l:ISpaans fc Silkll2006l: 
lOmukai et"al]|2008l: IWise et al.ll2008[ ). Indeed, the > 10"* K virial temperatures at the centers of lO'^ - 10^ M© halos heavily 
suppress H2 fractions, so it is a good approximation to take it that the opacity due to absorption of photons by H2 is also 
low. However, this does not imply that H^ fractions are negligible because it forms from reactants whose abundances are 
not suppressed by H2-dissociating backgrounds or high virial temperatures: 

R + e-^R-+-f, (A2) 

where 7 denotes the emission of a photon. In principle, the absorption of photons with energies > 0.75 eV by photode- 
tachment of H~ could be an important channel by which momentum can be imparted to inflow. To determine the rate 
at which momentum that could be transferred to the gas, it suffices to estimate the equilibrium rate of H~ formation, 
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since this is also the rate at which H will be destroyed and momentum will be acquired by the gas (see lAbel et al]|l997l : 
IChuzhov et al.|[2007l ). Using the density profiles from our calculations, which assume that the gas falls toward the star at 
the free-fall velocity above thii (Fig. 2), we have n(r) < 10^^ (r/rHii)~^ Numerical simulations of protogalactic 

collapse predict central free electron fractions /e ^ 10~^ (e. g. lShang et al.|[2"010l) . With these as upper limits along with 
the rate coefficient ka- for H~ formation, integration over the density profile of the halo yields the H~ formation rate 
Qh- outside the H ii region: 

Qh- < ka- I 47rr2/c^^^rfr ~ lO^^'s-^ , (A3) 



I-HII 



where fcn- ^ 10 cm -1 llWishartJIlQTllAbel et al.lll997^ . We assume that hydrogen is predominantly neutral outside 
the H II region and neglect the small H~ fractions that may form in the H ii region. Finally, to ensure that we have 
found a strong upper limit, we integrate this profile out to router — 10^^ cm, which is approximately the virial radius of 
the atomically-cooled halos at high redshift that host the supermassive stars we are studying. Even allowing for the high 
average photon energy for photodetachment /ij^h- ^ 10 eV for supermassive Pop III stars, the maximum momentum that 
could be transferred to the gas per unit time is 

^'^h-Qh- <1022gcms-^ (A4) 
c 

This is orders of magnitude smaller than the momentum of the accretion flow in equation 18, so destruction does not 
contribute to radiative feedback. 



Radiation Pressure from Trapped Lya and Balmer Series Lines 

In Section 2.2, we raised the possibility that Lya photons may slow down accretion by scattering from neutral atoms 
at the edge of the H ii region. While we showed that this is unlikely to alter the flow, mostly because Lya photons 
are confined to rLya, radiation pressure due to nebular emission lines in optically thick gas may be sufficient to reduce 
accretion. As also mentioned in Section 2.2, trapped Lya photons may be largely converted to Balmer series photons, 
in some cases after being destroyed by absorption by excited hydrogen or helium (jOsterbrock fc Ferlandl 120061 ). These 
Balmer photons can then propagate outward and efficiently cool the gas. The high temperatures, densities, and Lya 
photons trapped in the H ii region cause a large fraction of the hydrogen atoms to be excited to the n = 2 state there. 
From this state, either collisions or absorptions will further excite the atoms, which often later decay by emitting a Balmer 
series photon. 

To evaluate the degree to which radiation pressure from trapped line photons alters the dynamics of accretion, we first 
estimate the optical depth to Balmer photons both inside and outside the H ii region. Within rnii, the density of neutral 
hydrogen atoms can be estimated from equation 21 and assuming that n cx r"^, as implied by the nearly constant velocity 
profile of the gas in the H ii region (Fig. 2). The neutral hydrogen column density N-^ through the H ii region is then 

/•^'HII 

A^h(< fHii) ^ / fm{r)n{r)dr 

= /Hi(rHii) n(rHii)^ 2. 10^7 cm-2 , (A5) 

which, from the scaling in the second equation, can be shown to have only a weak dependence on the stellar mass for 
A'U = 10^ - 10^ Mq for our solutions to the minimum accretion rate. The cross sections for absorption of Balmer series 
photons by H in the n = 2 state are 19, 3.5, and 1.3 x 10~^^ cm^ for Ha, H/3, and H7, respectively. The higher energy 
Balmer series cross sections are all below these values. With these cross sections and N}i{< rnii), we can express the 
optical depth through the H 11 region solely as a function of the relative populations N2 / Ni of the n = 1 and n = 2 levels 
of neutral hydrogen in the H 11 region. We find that even for N2/N1 as high as ^ 0.05, the H 11 region is optically thin to 
all Balmer series photons, while the largest optical depth possible for Ha is thq — 20. However, this is still far below the 
optical depth for Lya photons, and Balmer series photons will eventually leak out of the H 11 region, although in some 
cases only after a number of scatterings. Therefore, while Lya photons in general will not escape the H 11 region, Balmer 
series photons will escape and allow the gas to radiatively cool, as we noted in Section 2.2. 

Next, we consider the optical depth to Balmer series lines outside the H 11 region. In this region, we set /hi = 1 and 
adopt a free-fall density profile n cx r^2 for the gas, as in Section 2.1. The neutral hydrogen column density beyond rnii 
is then 

/•oo 

A^h(> ran) ^ / n{r)dr 

' HII 

= 2n(rHii)rHii ^ lO'^ cm'^ , (A6) 

which again is not strongly dependent on stellar mass. 

Because the gas outside the H II region is well below the critical density at which excited levels in hydrogen would 
be populated and kept in equilibrium by collisions (and because Lya photons are mostly trapped in the H 11 region and 
cannot excite neutral hydrogen beyond run), the relative population of the n = 2 level of hydrogen is expected to be very 
small. In this case. iSchleicher et al., (,2010, ) find that in general N2/N1 < 10~^°. With this as an upper limit and the cross 
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sections for Balmer series photon absorption above, we find that the optical depths to Balmer series Unes are thq ~ 4 x 
10"^, th/j ^ 8 X 10"'^, and th^ ^ 3 x lO^'^. At such fow optical depths, Balmer series photons will propagate largely 
unimpeded through the inflow and, most likely, exit the halo with minimal effect on accretion. As these photons also move 
freely through the intergalactic medium (IGM), accreting supermassive stars probably have a distinctive observational 
signature, as we discussed in § 4. 

Thus, because outside the H ii region the optical depth to Balmer series photons is small and Lya photons cannot 
propagate across the H ii region boundary, we conclude that radiation pressure due to trapped emission lines will not 
be appreciable beyond run. However, it may be that this pressure is substantial within the H II region. To estimate its 
magnitude, we compare it to the ram pressure of the accretion flow. For the latter, we have 

-Pram = n^imuV^ 

dyncm"^ , (A7) 




100 km s 



-1 



where we have used equation 1 to express n in terms of v and Mace- We use the prescription of iBraun" fc Dekell (I1989D for 
the radiation pressure (see also lElitzur fc FerTand 1986): 

where Ni+i/ N\ is the ratio of the upper and lower level populations for the given transition, whose frequency is ftran- The 
line width Ai^tran 2 x 10^^ (T/10'* K)2 (In t)2 s^^, where r is the optical depth of the line. 

Assuming optical depths TLya = 10^ and TBaimcr = 20, both of which are rough upper limits in the H ii region, and 
Ni-^-il Ni = 1 in both cases, we find strong upper limits of 6 dyn cm~^ and 2 x 10~^ dyn cm~^ for the Lya and Balmer 
line radiation pressures, respectively. Comparing these pressures to Pram at thii, which is a lower limit for the H ii region 
due to its strong dependence on r in equation A7, we find that it is always at least a factor of ~ 4 larger than radiation 
pressure from optically thick lines. Therefore, lines will not play a large role in slowing accretion in the H II region. 

Gas Pressure 

Because infall is highly supersonic, gas pressure cannot decelerate the gas (jOmukai fc Inutsukall200^ . We have verified 
that this holds even when a gas pressure term is included in the equation of motion. We find that only for extremely 
high sound speeds corresponding to temperatures of at least 10^ K would gas pressure begin to impact the dynamics 
of the accretion f l ow. However, temperatures this high are not found in numerica l simulations of protogalactic collapse 
(jWise et al.ll2008t iRegan fc Haehneltll2009t IShang et al.]|2010t iJohnson et alll201lh or in H ii regions because the Balmer 
thermostat limits ionized gas temperatures to at most a few x 10''' K. 

Accretion Luminosity 

Our conclusions regarding the maximum mass of the supermassive star rest on the assumption that the critical accretion 
rate for a given stellar mass (equation 12) separates two regimes: below this rate, accretion is suppressed by radiative 
feedback and above this rate accretion proceeds. As shown in Section 2.1, it is clear that accretion is prevented for inflow 
rates below the critical rate because the ionization front breaks out to infinity and photoionization pressure halts infall 
at all radii. For accretion rates above the critical rate, the flow will arrive at the star with a velocity ^^(r*) > 0. This 
implies that the accreting material will have some kinetic energy that must be dissipated at the stellar surface, some 
portion of which will be radiation. This radiation could slow accretion within the H ii region by electron scattering or 
photoionizations . 

Let us assume that all of the kinetic energy of the flow is converted to radiation that propagates outward on impact 
with the star at some velocity v[r.^,). The Imninosity thus generated is 

L... = %^ , (A9) 

which yields a total momentum in photons of iacc/c = Macc^^^/2c. Comparing this to the momentum of the infalling 
gas, MiiccV, we see that the momentum in the accretion luminosity is a factor v/2c lower than that of the gas. Therefore, 
accretion radiation cannot halt the flow in the steady-state approximation. However, bouts of massive accretion could 
in principle generate enough radiation to slow infall at larg er radii where gas has a l ower momentum, like the episodic 
accretion onto black holes in the early universe discussed bv lMilosavlievic et al.l ( 20091 ). 



We note that an additional effect of accretion at very high rates (much hig her than those in equation 12) is the trappin g 
of radiation by e lectron scattering or oth er sources of opacity in the flow re.g. lBegelmanlll978l[2010l:IWvithe fc Loebll201l[ ). 
As discussed by lOmukai fc Fallal (j2003l see also Stabler et al. 1986) in the context of primordial protostellar accretion, 
radiation trapping can lead to rapid expansion of the stellar photosphere that halts the growth of the star. We note, 
however, that such expansion also causes the effective temperature of the star to fall with the increase in surface area of 
the photospher4I3i which in turn leads to a softening of the emitted radiation and a lower ionizing photon emission rate 
Q. This, and the fact that accretion rates greater than those in equation 12 result in smaller H II regions, suggests that 
the flow cannot be stopped, at least not by photoionization pressure. 

1 

At constant luminosity, the effective temperature at the stellar photosphere scales as T^g oc ' 



